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We propose new techniques to implement numerically the overlap-Dirac operator which exploit the physical 
properties of the underlying theory to avoid nested algorithms. We test these procedures in the two-dimensional 
Schwinger model and the results are very promising. We also present a detailed computation of the spectrum and 
chiral properties of the Schwinger Model in the overlap lattice formulation. 



The overlap ||l| formulation of lattice fermions 
provides a definition of a lattice Dirac operator D 
which avoids the doubling problem and preserves 
the relevant symmetries of the continuum theory, 
most notably chiral symmetry in the limit of van- 
ishing fermion mass. Unfortunately, this welcome 
development has come at a price: the numerical 
calculation of the matrix elements of the prop- 
agator D'^ and the inclusion of Det(D) in the 
measure entail a substantially increased computa- 
tional burden, which severely constrains the max- 
imum lattice size for viable simulations. In this 
talk we present results of an exploratory study 
for a numerical approximation which exploits the 
physical properties of the system. It proceeds 
through the projection over a subspace which 
has a substantially reduced number of degrees 
of freedom but still captures the relevant non- 
perturbative long range properties of the model. 

The Schwinger model is an ideal laboratory to 
test our approximation: it has many features in 
common with QCD, most notably chiral sym- 
metry, and it is simple enough that it is possi- 
ble to compute and Det(i3) exactly within 
the overlap formulation. Our primary goal is to 
compare the results of our approximation to ex- 
act results. Nevertheless the results of the exact 
calculations are interesting per se, because they 
are more extensive than what, to the best of our 
knowledge, has been obtained up to now and val- 
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idate in an impressive manner the advantages of 
Neuberger's formulation of lattice fermions. This 
talk is based on a work which will be the sub- 
ject of a forthcoming publication ||], where all the 
notations used here and the details of the simu- 
lations are reported. 

1. The Schvifinger Model with Overlap 
Fermions 

The Euclidean action of the Schwinger Model 
in the overlap regularization is given by 
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where U^v{x) is the standard Wilson plaquette, 
ij^i are Nf fermions fields with mass m and (5 = 
l/(ag)^, g being the bare couphng constant. D is 
the massless overlap-Dirac operator 
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where Dw is the Wilson-Dirac operator and a is 
the lattice spacing. The Neuberger-Dirac opera- 
tor satisfies the 75-Hermiticity condition 
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and the Ginsparg- Wilson (GW) relation 

-/5D-^ + D-^r. - a75 (4) 

The latter equation implies a continuous symme- 
try of the fermion action in Eq. at finite lattice 
spacing [| 
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where 75 = 75(1 — aD), which may be interpreted 
as a lattice form of the chiral symmetry. The 
chiral transformations in Eq. (j^) lead to Ward 
identities analogous to the continuum ones and 
the anomaly term arises from the non-invariance 
of the fermion integral measure 0|. The corre- 
sponding flavor non-singlet chiral transformations 
are defined including a flavor group generator in 
Eq. (H). They imply the Axial Ward identities 
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where the axial current and the pseudoscalar den- 
sities are 
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with tp = {jpi, . . .tpff^.), and A'^ are the generators 
of the SU{Nf) flavor group. In two dimensions 
the Axial current is related to the Vector one as 
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where e^^ is the totally antisymmetric tensor. 

In the Nf — 2 massless model, there is a triplet 
of massless particles (tt) and one massive singlet 
particle (77) with mass (M^/g)^ = 2/-K. The 
correlation functions of the vector currents are 
the same as for free particles. For small masses 
the corrections to the massless results can be ob- 
tained by "chiral perturbation" theory |^ which 
to first order gives 
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where 7 — 0.577216 ... is Euler's constant. 
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Figure 1. M^^j g vs (mj g)^^"^ for the full operator 
for Nj = 2. 



2. Numerical results vifith the full operator 

In order to test our approximation method and, 
at the same time, to increase the body of in- 
formations on the lattice Schwinger model, we 
performed an extensive simulation with the ex- 
act Neuberger-Dirac operator. We considered the 
Schwinger model at /3 = 6 on a 24 x 24 lattice. On 
a lattice of this size, the discretized Dirac opera- 
tor is a complex matrix of dimension 1152 x 1152, 
for which we could use full matrix algebra subrou- 
tines without excessive burden on the resources 
available to us (Boston University SGI/Cray Ori- 
gin 2000 supercomputer). The masses have been 
selected on the basis of previous results which in- 
dicated that, for most of the values we were plan- 
ning to consider, the lattice would span at least 
a few correlation lengths. We generated 500 in- 
dependent configurations of the gauge variables 
Ufi{x) distributed according to the pure gauge 
Wilson action of Eq. (^. For calculations with 
one and two flavors of dynamical fermions, we in- 
corporated the determinant of the lattice Dirac 
operator in the averages of the correlation func- 
tions 1^ . While with large variations of the deter- 
minant this way of proceeding would lead to an 
unacceptable variance, in the present calculation 
we found the range of values taken by Det(D) to 
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Figure 2. m^^^ jg vs mjg for the full operator 
for Nf = 2. 



be sufficiently limited to warrant our averaging 
procedure (see Fig. ^). 

For each configuration, we performed a singular 
value decomposition 
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where C/, JJ are unitary matrices and A is diag- 
onal, real and non-negative. The operator V = 
X/y/lOX in Eq. | is then given by 
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We proceeded then to the diagonahzation of V , 
calculating all its eigenvalues and eigenvectors. 
From these, it is straightforward to calculate both 
the determinant of the Neuberger-Dirac operator 
D as well as its associated propagator for any 
value of the fermion mass m. The full diagonal- 
ization of V is computationally more demanding 
than the direct calculation of Z?^^, which typi- 
cally gives also Det(D) as a by-product, but we 
were interested in the actual spectrum of V for 
comparison with the approximations that will be 
discussed later. From the fermion propagators we 
calculated the meson propagators projected over 
zero momentum. We focused on the two point 
correlation functions of the singlet vector current 
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Mesons and quark masses for Nf — 2. 





f3 = 6.0, V 


= 24^ Nf = 2 




m/g 




m-^/g 


m^/g 


Analytic 











0.7979 


Full Operator 


0.00 


-0.004(2) 


-0.001(65) 


1.00(25) 


0.0244 


0.0233(2) 


0.14(3) 


1.0(5) 


0.0485 


0.0468(6) 


0.26(5) 


1.0(4) 


0.0960 


0.0945(13) 


0.44(4) 


1.1(2) 


0.1427 


0.1420(14) 


0.57(3) 


1.1(2) 


0.1884 


0.1890(13) 


0.69(2) 


1.2(1) 


0.2333 


0.2353(11) 


0.79(1) 


1.23(9) 


Fourier Approximation 


0.00 


-0.002(2) 


-0.005(66) 




0.0244 


0.018(1) 


0.11(3) 




0.0485 


0.044(1) 


0.24(5) 




0.0960 


0.092(1) 


0.43(4) 




0.1427 


0.139(1) 


0.57(3) 




0.1884 


0.184(1) 


0.68(2) 




0.2333 


0.228(1) 


0.79(1) 




Multi-grid Determinant 


0.00 


-0.004(1) 


-0.002(40) 


1.09(25) 


0.0244 


0.0234(1) 


0.14(2) 


1.1(4) 


0.0485 


0.0469(4) 


0.26(3) 


1.1(4) 


0.0960 


0.0946(9) 


0.44(2) 


1.1(2) 


0.1427 


0.1420(10) 


0.57(2) 


1.2(2) 


0.1884 


0.1889(9) 


0.68(1) 


1.2(1) 


0.2333 


0.2351(8) 


0.79(1) 


1.26(9) 



and the triplet current 
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because they are saturated by a single particle 
contribution in the massless limit. We added to 
the averages the correlators obtained from the in- 
terchange of X and y in Eq. |l^ for a further gain 
in statistics. 

From fits to the meson correlators we extracted 
the meson masses in a standard manner. Our 
results for the meson and quark masses as func- 
tions of the fermion mass parameter of the ac- 
tion arc reported in Figs. |l| and ^. The values we 
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Figure 3. Structure of the Wilson operator in 
momentum space and in a smooth gauge. 



obtained for the masses are also reported in the 
Table. The errors have been evaluated with the 
jackknife method. Fig. ^ illustrates the behav- 
ior of the isotriplet mass A/^ in the model with 
two flavors. This mass is expected to vanish for 
m = and chiral perturbation theory predicts an 
m^/'^ dependence on m (see Eq. ^. To test this 
functional dependence we performed two different 
fits. First the pion masses were fitted as 

2/3 
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using only the highest four quark masses for 
which 24 X Mt, > 4. The results of the fit 
A = -0.001(65) and B = 2.10(14) confirm that 
Neuberger's fermions preserve chiral symmetry 
and give values of B in very good agreement with 
the theoretical prediction. The dashed line in 
Fig. |l| is the result of the fit. We then fitted the 
pion mass as 

^=B('l^V (14) 
9 \9 J 

obtaining B = 2.10(17) and 6 = 0.67(6) in per- 
fect agreement with the expected functional de- 
pendence. Equally gratifying is the comparison of 



the value of the fermion mass m in the Lagrangian 
with the value m^^^ that can be extracted from 
the axial Ward identity Eq. ^ 

In the Table we also report results for the sin- 
glet mass Mrj versus m'*/^, always for Nf = 2. 
The singlet meson propagator is given by the dif- 
ference of the connected and disconnected terms 
in the correlator and, because of the cancellations, 
the errors are larger than in the triplet case. The 
numerical results are, however, consistent with 
the theoretical prediction of Eq. ||. 

3. Approximation in the Fourier Space 

The long range physical properties of the 
Schwinger Model should be determined by the 
eigenstates of D with eigenvalues in the neigh- 
borhood of A = m (i.e. the eigenstates of V with 
eigenvalue Xy closest to —1), since these are the 
states with are expected to go over the physical 
states of the continuum in the limit a —^ 0. This 
suggests that it should be possible to reconstruct 
the physical observables from these states alone. 
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Figure 4. Comparison of the V2-triplet corre- 
lation function computed with the full operator 
and the Fourier approximation on a configuration 
with (3 = 0. 
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Figure 5. Comparison of the V2 -triplet corre- 
lation function computed with the full operator 
and the Fourier approximation on a configuration 
with (3 = 4. 
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Figure 6. M^^/g vs [m/gY/^ for the full oper- 
ator (circles, blue) and the Fourier approximation 
(stars, red) for Nf = 2. 



Here we would like to make the point that the 
overlap formulation is particularly well suited for 
the implementation of the above approximation, 
since the unitarity of V provides, in some sense, 
a constraint on the form of the approximation 
itself. Of course, one must be careful in attempt- 
ing any approximation based on neglecting the 
short- wavelength part of the spectrum, even if the 
corresponding states are largely lattice artifacts, 
since one knows that in quantum field theory the 
infrared and ultraviolet components of the spec- 
trum are subtly related. Thus, the chiral eigen- 
states with Ay = — 1 ("zero modes"), which V 
exhibits in presence of a gauge field with non- 
trivial topology, find their counterpart in states 
with opposite chirality and Ay — 1. 

Here too, however, the special features of the 
overlap formulation come to the rescue, since, as 
we will show, both the presence of zero modes and 
the correspondence between Ay = — 1 and Ay = 1 
eigenstates of opposite chirality is preserved, as 
we found in our study, by the projection over a 
subset of physical states. 

One possible scheme of approximation which 
is computationally very convenient consists of 



performing a projection over states of low mo- 
mentum in Fourier space, after gauge fixing 
to a smooth gauge field configuration^. In a 
smooth gauge, because of the suppression of 
short-wavelength fluctuations due to asymptotic 
freedom, one expects the structure of the Wilson 
operator to become more and more diagonal in 
momentum space for increasing momenta. This 
notion is schematically illustrated in Fig. ^ and 
also underlies the technique of Fourier accelera- 
tion for the calculation of quark propagators . 
Accordingly, we implemented the following ap- 
proximation. For each configuration we factored 
out the topology, if any, as described in ||] and 
fixed the Landau gauge on the resulting configu- 
ration by demanding that the function 

G = J2 MU^^ix) + U^{x - afi)] (15) 

be maximal. A relaxation procedure produces 
several local maxima (Gribov copies). Among all 
these configurations, we selected the one that pro- 
duced the maximum for G, subject to the further 
constraint that for all x and /i JieU^{x) > 0.5. 

similar but gauge invariant procedure is described in 
which, however, turns out to be numerically more ex- 
pensive. 
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Figure 7. Sketch of the blocking procedure. 
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Figure 8. Det(D) on the blocked lattice versus 
Det(Z3) on the full lattice. 



Eventually the gauge transformation obtained 
was applied to the original configuration. 

Once the gauge has been fixed, the vector 
space on which the fermion matrix D acts is di- 
vided into two parts: a long range (LR) sector 
spanned by the vectors of the Fourier basis \pi) 
(i — 1 , . . . n) with pf < p\ and its complemen- 
tary short range (SR) subspace. In the LR sector 
the matrix elements of the Wilson operator are 
computed and Neuberger's operator is obtained 
with the unitarity projection in Eq. and in- 
verted. In the SR subspace the propagator is ap- 
proximated by its free form. The approximated 
operator is 75-Hermitian and it satisfies the GW 
relation. 

In Figs. U and (||) we compare the triplet meson 
correlation functions Ct (t) obtained with the full 
Neuberger operator and with the approximated 
one for a configuration with trivial topology and 
for one with the topological charge Q — 4. In 
Fig. the triplet meson masses for the full solu- 
tion and the Fourier approximation are reported. 
The LR subspace used is 1/4 of the full space. It 
is interesting to note that the approximated so- 
lution can be improved if used as preconditioning 
of the standard algorithms. 



4. Multi-Grid Approximation for the De- 
terminant 

In the spirit of the long range approximation, 
another possibility that can be explored is the 
coarse graining of the gauge field 1^. Start- 
ing from an arbitrary gauge configuration on an 
N X N lattice, we apply the blocking procedure 
sketched in Fig. ^to obtain a gauge configuration 
on an Y X Y If^ttice, which carries all the relevant 
long range features of the large lattice. 

The original lattice is divided into ^x^- blocks 
of 2^ sites. Each block, after the local Landau 
gauge is imposed, corresponds to a site of the 
coarse lattice. The remaining links are averaged 
as in Fig. ^ to give the corresponding links of the 
coarse lattice. 

The absence of additive mass renormalization 
simplifies the definition of the Neuberger operator 
on the coarse lattice because no fine tuning of 
the mass parameter is needed. We calculated the 
fermion determinant on the blocked lattice and 
compared it to the one of the full lattice as shown 
in Fig. ||. 

We incorporated the blocked determinant in 
the averages of the correlation functions. We 
compare these results with those obtained with 
the determinant of the full operator. For the pion 
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masses the agreement is remarkable (see Fig. ||). 
We did the same for the massive Schwinger boson 
and the results are reported in the Table. Also in 
this case, even if the observables are more noisy, 
the agreement between the results is remarkable. 

The appHcation of these ideas to speed up 
QCD simulations with Neuberger fcrmions is in 
progress. 
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